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We measured the angular dependence of the c-axis magnetoresistance p c (B) of Bi2Sr2CaCu208+a 
irradiated with heavy ions. At temperatures near 68 K the scaling of p c (B) with the c-axis magnetic 
field component B± is broken and the in-plane field, Bu, affects p c . At this temperature, at a 
specific field B cr « B&/2, magnetoresistance becomes independent of field orientation. This crossing 
point allows us to estimate the correlation length L of pancake positions along the c axis. We 
find L/s ~ fOO at B = B cr , where s is the interlayer spacing. This provides evidence of strong 
enhancement of pancake alignment in the vortex liquid in crystals with columnar defects. 
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The properties of the vortex liquid in high-temperature 
superconductors, HTS, in the presence of strong disorder 
is one of the most challenging problems in the physics of 
the vortex state. It is well established now that in pris- 
tine Bi 2 Sr2CaCu208+(5, Bi-2212, crystals vortices form 
a weakly c-axis correlated pancake liquid. Thermody- 
namics and intralayer dynamics of this liquid are defined 
mostly by pancake concentration, i.e., by the c-axis mag- 
netic field B± 0. This behavior is however, strongly af- 
fected by the presence of correlated disorder. The most 
effective pinning centers are produced by heavy ion ir- 
radiation. Irradiation of the HTS by energetic ions pro- 
duces amorphous tracks, where the superconductivity is 
suppressed. Such columnar defects, CDs, with radii com- 
parable with the superconducting coherence length, are 
ideal for pinning vortex lines, but their effect on pancake 
vortices is less obvious. The main questions are: a) are 
pancake vortices positioned mainly inside columnar de- 
fects in the liquid phase?, and b) are c-axis correlations 
enhanced in the presence of CDs or do they remain sim- 
ilar to that in the liquid phase of pristine crystals? 

Important information about the effect of CDs on the 
pancake liquid was obtained from reversible magneti- 
zation measurements [^j. In pristine crystals reversible 
magnetization M monotonically increases with the mag- 
netic field B±. In irradiated crystals, due to gain in pin- 
ning energy, penetration of vortices into crystals becomes 
more favorable than in pristine crystals, and the diamag- 
netic moment drops in the presence of CDs at low fields 
B± <C -B$ when all vortices may occupy CDs. However, 
as B± increases towards the matching field intersti- 
tial vortices start to appear, losing the advantage in pin- 
ning energy associated with the CDs. This results in a 
drop of magnetization (i.e., in an increase in the diamag- 
netic moment) in the field interval between rs B$ /4 and 
B$. At larger fields the difference between pristine and 
irradiated crystals almost vanishes. Such an anomaly in 



M(B±) dependence was observed in the vortex liquid up 
to rather high temperature, indicating that mobile pan- 
cakes, inherent to the vortex liquid, are localized largely 
onto CDs, even at temperatures close to T c . Recent study 
of the vortex-lattice melting in weakly irradiated Bi-2212 
also confirmed this picture ||. However, the thermo- 
dynamics of the vortex state depends weakly on c-axis 
correlation of pancakes; it is determined mainly by in- 
tralayer interactions, which are much stronger that those 
associated with the c-axis vortex structure |[|,|| . For this 
reason reversible magnetization is determined by pancake 
concentration and scales with B±. Thus magnetization 
measurements show that pancakes are positioned mainly 
inside CDs at B± < B$/3 but they do not provide infor- 
mation on c-axis correlations of pancakes inside CDs. 

In contrast, Josephson interlayer properties of highly 
anisotropic Bi-2212 are extremely sensitive to the c-axis 
correlations of the pancakes, because Josephson current 
depends on the interlayer phase difference. Pancakes, 
aligned along the c-axis do not contribute to the phase 
difference, but those, shifted due to thermal fluctuations 
or pinning do. This leads to larger phase difference in 
the uncorrelated liquid and thus to the suppression of 
the c-axis superconducting current. Measurements of 
Josephson plasma resonance, JPR, reveal for the first 
time that c-axis correlation in the vortex liquid in the 
presence of CDs depends on B± nonmonotonically in the 
temperature interval 60 - 69 K, showing enhancement 
at B± ~ _B$/3 0], while in the liquid phase in pristine 
crystals correlation drops with B± at all temperatures. 
Nonmonotonic behavior with B± was found also in the 
dependence of the c-axis resistivity, p c , which is also de- 
termined by interlayer Josephson current Q. Namely, 
p c in the same temperature interval exhibits an increase 
with B± at B± <C -S$, a flattening or even a decrease 
in the field interval between 0.2 B$, and 0.4 B$, followed 
by further increase of p c at higher B±. Unfortunately, 
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neither JPR nor p c (B±) measurements provide sufficient 
information to estimate the enhancement of the c-axis 
correlation quantitatively. 

It was shown in |||| that dependence of JPR frequency 
and p c on the parallel component of the magnetic field 
provides complete information about c-axis phase corre- 
lations. Based on this idea, in this paper we study the 
angular dependence of p c at different orientation of the 
magnetic field B in the correlated pancake liquid phase 
in the presence of CDs. We show that scaling with B± 
is absent in this phase. Moreover, we found, that at 
B = B cr k, B&/2, this angular dependence is accurately 
compensated, providing a crossing point where p c is in- 
dependent of the orientation 9 of the magnetic field with 
respect to the crystal c-axis over a wide range of 9. Such 
an anomaly gives an opportunity to estimate the corre- 
lation length of pancake positions along the c axis at the 
crossing field. 

For our experiments high quality Bi-2212 crystals 
(T c ~ 85 K) of about 1 x 1.5 x 0.02 mm 3 were used. The 
irradiation by 1.2 GeV U 238 -ions was performed on the 
ATLAS accelerator (Argonne National Lab.). According 
to TRIM calculations these high energy ions produce in 
Bi-2212 crystals continuous amorphous tracks with diam- 
eter 4-8 nm and length 25-30 /im. Below we present the 
results for samples irradiated with an effective density 
of columnar defects corresponding to the matching field 
B$ = 1 T. Another pristine sample was used as a refer- 
ence. We checked also several crystals from a different 
synthesis for universality of the obtained results. 
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FIG. 1. Dependence of c-axis magnetoresistance p c on 
c-axis field for two crystals: nonirradiated, B$ = 0, and ir- 
radiated with effective density of CDs B$ = IT. The latter 
displays dip due to enhanced vortex correlations in the pan- 
cake liquid. Insert: schematic experimental set-up. 

Our measurements were performed using a warm bore 
insert in a liquid helium cryostat with a 9 T supercon- 
ducting magnet, and also in a liquid nitrogen, LN, dewar 
installed into an electromagnet that provides a magnetic 
field up to IT. The sample was attached to a sample 
holder with goniometric stage, allowing rotation of the 



sample with respect to the applied field with an accu- 
racy of about ±0.1°. The sample was immersed in a LN 
bath in order to provide accurate temperature stabiliza- 
tion and to avoid heating effects. Magnetic fields applied 
along the c-axis of the sample were monitored by a Hall 
sensor attached to the sample stage. Two pair of silver 
contacts pads were deposited on both sides of the sam- 
ples for standard 4-probe transport measurements and 
fine gold leads were attached using silver epoxy. The re- 
sistivity for the contact pair at room temperature was 
about 1 — 3 fi. We measured R c as a function of applied 
magnetic field B at different angles 9 with respect to the 
c-axis of the crystal (see insert in Fig. 1). 

In Fig. 1 we present the normalized magnetoresistance 
Pc = Pel Pn as a function of applied field along the c- 
axis measured at T — 68 K. Here p n is the normal state 
resistivity at T — 120 K. After the sharp onset at the 
irreversibility line, p~ c for irradiated crystals displays a 
pronounced dip due to enhanced vortex correlation asso- 
ciated with filling of the CDs reported in Ref. 0. At 
fields B ^> Bg, the magnetoresistance approaches the 
normal state value. In the pristine crystal the resistance 
increases with B± monotonically. 
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FIG. 2. Dependence of p c versus applied field B for differ- 
ent 8, progressively increased with step 5° (left). The same 
data plotted as a function of c-axis field B c (right). The scal- 
ing of p c with B± is broken in the correlated vortex liquid. 

In the left panel of Fig. 2 the family p c vs. B at differ- 
ent angles is presented. The magnetoresistance decreases 
with the angle at B > _B$/2 (as in the pristine crystals), 
reflecting decrease of B±, but increases with the angle 
at low fields (0.1B$ < B < B&/2), in the region where 
reentrant enhancement of c-axis correlation was observed 
in Pc(B±) dependence. 

In our experiment the inclination of the sample in the 
applied field results in two effects: decrease of the c- 
axis field B± = B cos9 and increase of the in-plane field 
£?H = B sin9. In order to separate these two effects we 
present p c as a function of the c-axis component B±_ in 
the right panel of Fig. 2. The magnetoresistance of the 
sample scales with B± at low and at high fields. However, 
in the range B&/A < B < this scaling is broken. 
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This anomalous behavior occurs in the range of field and 
temperature where enhanced vortex correlations in the 
pancake liquid were observed in JPR and in p c (B±). 

In Fig. 3 we blow up the region where the scaling is 
broken. It is clearly seen on the right panel that p~ c here 
depends on the in-plane field, increasing gradually with 
Bu . The most outstanding feature can be seen on the 
left panel of Fig. 3. Here p c curves at different angles are 
presented as a function of B. They all cross at a single 
point, which we denote as B cr . This crossing point is ob- 
served for all the crystals which display the dip of p c . It 
means, that at B = B cr , angular dependence of magne- 
toresistance practically vanishes, i.e., decrease of p c due 
to decrease of B± is accurately compensated by increase 
of the magnetoresistance with increasing Bu. This ef- 
fect can be explained in the framework of the approach 
developed in ||[| and will be used below to estimate 
quantitatively the characteristic length Ri of the phase 
difference in-plane correlation function. 




FIG. 3. Blow up of the anomaly for the irradiated crys- 
tal. Dependence of p c versus applied field B (left) and versus 
c-axis field B± (right). 6 increases from to 45° with step 
5°. The crossing point B cr corresponding to the absence of 
the angular dependence is marked on the p c vs. B panel. 



Let us first discuss how the in-plane field component 
affects the c-axis magnetoresistance. An approach de- 
veloped initially for the field behavior of JPR frequency 
was extended recently for c-axis transport properties 
In Josephson coupled superconductors in the pres- 
ence of a c-axis current, the voltage V n , n +i is induced by 
slips of the phase difference <pn,n+i{r, t) between the lay- 
ers n and n + 1, as described by the Josephson relation 
V n>n +i — (h/2e)(p ntn +i. Here r = x,y are coordinates in 
the ab plane, and t denotes the time. The c-axis conduc- 
tivity in the vortex liquid state er c = l/p c is determined 
by the Kubo formula |l0| 
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<T c (B ± ,B ll ) = (sJ*/2T) dtJdrS(r,t), (1) 

S(r,t) = 2(sin^ n , n+ i(0,0)sin<p„ jn+ i(r,t)) 

w (cos[^„, n+ i(r,i) - <p n;n+1 (0,0)]), (2) 



where j7o is the Josephson critical current, s is the inter- 
layer distance, and < . . . > means thermal average and 
average over disorder. Time variations of the phase differ- 
ence are caused mainly by mobile pancakes | JlT[ | induced 
by -Bj_, while the parallel field component induces a sta- 
tionary phase difference in the lowest order in Josephson 
coupling. We split (r, t) - i/> n ,„ + i(0, 0)] into the 

contribution induced by pancakes and that caused by the 
unscreened parallel component B\\ . Assuming that Bu is 
along the x axis, we obtain: 



9V„+i(0,0) - y> n ,„+i (r,t) = 

[Vn,n+l(0,0) - cpn tn+1 (r,t)] Bu 



(3) 
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In a single point-like junction the phase difference in- 
duced by B\\ results in the Fraunhofer pattern of Joseph- 
son critical current as a function of the magnetic field 
parallel to the junction. In our case of a multilayer super- 
conductor the phase difference induced by Bu interferes 
with that induced by pancakes. Thus we obtain 



a c (B ± ,Bu) = (nsji/T) / drrG(r,B ± )J (aB«r), (4) 



G(r,B ± )= / dtS{r,t 7 B ± ), 
Jo 

where Jq(x) is the Bessel function, a = 2irs/&o ($o is 
a flux quantum), and the function G(r,B±) depends on 
correlations of the phase difference induced by pancake 
vortices. For small Bu (small angles) we expand the 
Bessel function in Bu : 
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<r (Bj.,B||) wcr c (Bj.,0) 
Rf(B±) = [ drr 3 G(r,B ± )/ I drrG(r,B±) 



(5) 



Here the in-plane correlation length Ri(B±) describes 
the decay of the phase difference correlation function. 
When pancakes are positioned mainly inside CDs, the 
characteristic length of this decay, R\, gives direct in- 
formation on the c-axis correlation of pancake positions 
because drop of the phase difference correlations in the 
ab plane is caused by interruptions in the pancake ar- 
rangement along CDs ||. The characteristic length, L, 
of the pancake density correlation function is related to 
i?i as L/s « Rl/lOa 2 , where a = ($ /-B ± ) 1/2 is the in- 
tervortex distance || . This expression and Eq. (||) are a 
key points for further discussion. 

It is clear from Eq. (||) that generally a c depends upon 
both components of the field. However, in an uncor- 
related liquid, when i?i w a, the effect of the in-plane 
field is small and can be observed only in high fields 
Bu $o/sa. As a result, in the uncorrelated pancake 
liquid the c-axis conductivity scales with B± in fields 
B <C <f>o/ sa - For an irradiated sample, as filling of CDs 
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progresses, vortices start to form stacks and the correla- 
tion length Ri significantly exceeds a, reaching a max- 
imum value near B$/3. Here the effect of the in-plane 
field becomes significant and scaling of magnetoresistance 
with B± breaks. As the field B± further increases, the 
fraction of the interstitial vortices increases. These vor- 
tices introduce additional disorder to the system which 
results in decay of correlations, in reduction of R\ down 
to ~ a and, consequently, the effect of Bu drops. Then, at 
elevated fields above Bg, scaling of p c with B± is restored. 
This scenario describes qualitatively well our experimen- 
tal results. 

Now we show that the crossing point in Fig 3 allows us 
to estimate the correlation radius R\(B cr ). Using Eq. (||) 
for the variation of a c at small angles we expand as: 



8a c (B h B ± ) 



d 2 a c (5B n ) 2 da c 



dB 2 



dBj 



5B< 



(6) 



Substituting Bn rs BO and B±_ ~ 
5a c (B h B ± )*lB(^B 



B(l - 9 2 /2) we obtain 
da, ' 



dB, 



(sey 



(7) 



To conclude, we have presented evidence for the pres- 
ence of a partially aligned vortex liquid in irradiated Bi- 
2212 from the angular dependence of c-axis magnetoresis- 
tance. In the range of temperatures and magnetic fields 
where c-axis correlations develop, interlayer transport be- 
comes much more sensitive to the in-plane component 
of the magnetic field in comparison with the uncorre- 
cted liquid in pristine crystals. Then scaling of p c with 
B± breaks, and, at the field B cr , angular dependence 
becomes very weak. From this crossing point we esti- 
mate the c-axis correlation length of pancake positions 
L/s w 100 at B± w £?$/2. The origin of this alignment 
can be due to magnetic intralayer interaction of pancakes 
which favors similar filling of CDs due to their geometry. 
However an additional effect of interlayer magnetic inter- 
action of pancakes cannot be excluded so far. 
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Independence of a c on the angle 9 at B cr occurs when 
the expression in the brackets becomes zero. With the 
help of Eq. (H) we obtain the correlation radius R\ : 



R 2 (B C 
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2B C 



8\na c (B ± ,0) 



dB\ 



(8) 
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From the data presented in Fig. 3, we obtain 
(lncr c (S_L,0))' Ri -23 T" 1 at B = B cr w 0.47 T. Using 
Eq. (^) we calculate R\/a Ri 32. This value is in rea- 
sonable agreement with the result obtained in [^j , where 
Ri ~ 10a was found. Thus the correlations length L of 
pancakes positions along the c axis is Ri 100s. 

We would like to emphasize the difference between our 
observation and results obtained in flux-transformer ge- 
ometry p2|| . In those latter experiments current was ap- 
plied to the top surface of the sample and voltage was 
measured in the both top and bottom layers. In some 
range of the magnetic fields and temperatures top and 
bottom voltages coincide, indicating a similar motion of 
vortices in all layers. This was considered as evidence 
of pancake coupling along the c-axis. Unfortunately, the 
situation here is rather complicated because of mixing 
of the in-plane and c-axis resistivities. Namely, the cur- 
rent applied along the top layer penetrates deep into the 
sample along the c axis and flows along the bottom layer 
as well as along the top layer. This results in a similar 
motion of the pancakes through the full sample thick- 
ness, even without c-axis correlation provided p c is small 
enough. As we have shown, the c-axis transport measure- 
ments as a function of in-plane field component are free 
of this additional effect and are very sensitive to pancake 
correlation between adjacent layers. 
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